Transport and percolation in a low-density high-mobility two-dimensional hole system 
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We present a study of the temperature and density dependence of the resistivity of an extremely 
high quality two-dimensional hole system grown on the (100) surface of GaAs. For high densities 
in the metallic regime (p ^ 4 x 10® cm~^), the nonmonotonic temperature dependence (~ 50 — 300 
mK) of the resistivity is consistent with temperature dependent screening of residual impurities. At 
a fixed temperature of T= 50 mK, the conductivity vs. density data indicates an inhomogeneity 
driven percolation-type transition to an insulating state at a critical density of 3.8 x 10® cm~^. 



The low temperature conduction properties of an in- 
teracting two-dimensional (2D) electron gas in the pres- 
ence of disorder remains a topic of fundamental interest 
and some controversy. It has been understood for some 
time that all non-interacting electronic states are local- 
ized in the presence of disorder at r=0 [ij. This model 
was challenged by the observation of Kravchcnko and co- 
workers 01 of an apparent mctal-to-insulator transition 
in Si MOSFETs. Subsequently a variety of models have 
been put forth 0, H, 0] to explain a wealth of experi- 
mental data. The 2D metal-to-insulator transition (MIT) 
has now been studied in a variety of semiconductor sys- 
tems 13, 0, i, i, M, El, El, [13, El- Yet the question 
remains: does the transition from insulating to metal- 
lic behavior observed in certain 2D systems represent a 
quantum phase transition to a novel metallic state sta- 
bilized by interactions or can the experimental data 
in the metallic regime be explained in terms of finitc- 
temperaturc Fermi liquid behavior [3]? The MIT in Si 
MOSFETs has been interpreted as evidence for quantum 
critical behavior Q. However, recent experimental work 
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15j in ultra-high mobility n-GaAs 2D systems has 
interpreted the MIT as a density inhomogeneity induced 
percolation transition. 

In this Letter, we present a study of the transport prop- 
erties of a clean two-dimensional hole system (2DHS) in 
the vicinity of the putative metal-to-insulator transition. 
Our carbon-doped 2DIISs display mobilities in excess of 
2x10^ cm^/Vs at densities p<3xl0^*' cm^^, the highest 
ever achieved in a 2DHS in this density range. The use of 
extremely high mobility samples allows us to probe the 
regime of large {rs~Ec /Ep, where Ep is the Coulomb 
interaction scale and E^? is the Fermi energy) in a clean 
system, where the effects of carrier-carrier interactions 
should be most apparent. Interestingly, we find that 
for high densities in the metallic regime the nonmono- 
tonic temperature dependence of the resistivity is con- 
sistent with temperature dependent screening of resid- 
ual impurities. At a fixed temperature of T=50 mK, 
the conductivity vs. density data indicates a disorder- 
induced percolation-type transition to an insulating state 
at low density. While the critical density (pc ~ 3.8 x 10^ 



cm~^) and the interaction parameter (r^ ~ 45) for our 
sample are among the lowest and highest ever respec- 
tively recorded in the 2D MIT literature, we find that 
the data are well explained within this simple frame- 
work. The fact that our critical density is similar to 
the critical density for 2D MIT found in ref. sug- 



gests that the effective 2D MIT in GaAs electron and 
hole systems at experimentally accessible temperatures 
is a disorder-driven crossover phenomenon and not an 
interaction driven quantum phase transition, since the 
Ts value for our hole system is four times larger than the 



corresponding Vs in the electron system of ref. [11 1 



Our experiments utilize recently developed structures 
in which the 2DHS resides in a 20nm GaAs/AlGaAs 
quantum well grown on the (100) surface of GaAs that 
has been modulation-doped with carbon [l(i |. A back- 
gated Hall bar geometry is used in this experiment to 
tune the density from 2.9x10^° cm~^ to 2.9x10^ cm~^ 
in single sample. Cyclotron resonance studies of simi- 
larly grown carbon-doped quantum wells indicate that 
the hole effective mass is large, to^ ~ O.Stoc, where rUe 
is the free electron mass [TtI . To give a comparison with 
other 2D systems studied in the literature, we mention 
that the 2D n-GaAs samples of refs. 13, III have an elec- 
tron mobility of 5 x 10^ cm^/Vs at the carrier density of 
3 X 10^° cm~^ compared with our hole mobility of 2 x 10^ 
cm^/Vs at the same density. Taking into account the ef- 
fective mass difference (m = 0.07me for electrons in refs. 
[lo| and [m and m w 0.3me for holes in our samples) 
and remembering that mobility ^ cx r/m, where r is the 
transport relaxation time, we conclude that our 2D hole 
samples are comparable in quality to the corresponding 
2D electron samples of refs. [llj while is signif- 
icantly larger. In this manner the data presented here 
probe a regime inaccessible to the experiments of refs. 

Figure 1 displays magnetotransport at r=:50 mK. The 
exceptional quality of these samples is evident in the 
strength of the fractional quantum Hall state at iy=2/3 at 
B'^1 Tesla and the nascent features observed at v=5/3 
and 4/3 at even lower values of magnetic field. Figure 
2 displays the temperature dependence of the resistiv- 
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FIG. 1: Longitudinal and Hall resistance at zero gate bias. 
The fractional quantum Hall state at i/=2/3 is fully developed 
at B~1T while weaker features at j^=5/3 and 4/3 are also 
visible. 



ity at B=0 for several values of density p<9xl0^ cm~^. 
At all but the lowest densities, the temperature depen- 
dence is strongly nonmonotonic [3. Starting from the 
high temperature regime (T'^l K) the resistivity first in- 
creases vifith decreasing temperature, reaches peak value, 
and then decreases rapidly with further reduction in tem- 
perature. At p=9.0x 10^ cm~^, the low temperature drop 
in resistance is a factor of 2. As the density is reduced 
this crossover point, where dp/dT changes sign, moves 
to lower temperature, presumably tracking the reduction 
in the Fermi energy of our 2D system. The metallic be- 
havior persists to very low density, at p=4.1xl0^ cm~^ 
dp/dT still changes sign near T=50 mK. At p=2.9xl0^ 
cm~^, the lowest density accessed in our study, p vs. T 
is monotonic, suggesting that the sample is approaching 



an insulating regime. Yet at p=2.9xl0 cm~ the con- 
ductivity is not exponentially suppressed with decreas- 
ing temperature as expected for an insulating state, but 
rather, the conductivity displays a linear dependence on 
temperature (i.e. p oc T"^) in the T = 50—450 mK range 
[l9| . In this regime of low density, the Fermi temperature 
is as low as Tp ~ 250 mK. The impact of carrier degener- 
acy and the temperature dependence of ionized impurity 
scattering must be considered as discussed below. 

In order to understand the low temperature 2D trans- 
port properties in the metallic regime we have carried out 
a microscopic transport calculation using the Boltzmann 
theory In our calculation we assume that the 2D 
carrier resistivity is entirely limited by screened impurity 
scattering. We neglect all phonon scattering effects be- 
cause our theoretical estimate shows phonon scattering 
to be negligible for 2D holes in GaAs structures in the 
T <\K regime of interest to us. Formally, this approach 
is well-controlled in the regime of small Vg and kpl ^ 1, 
where kp is the Fermi wavevector and / is the mean free 
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FIG. 2: Temperature dependence of the resistivity for 50 
mK<T<l K for densities ranging from 9.0x10^ cm^'^ to 
2.9x10^ cm"^ 



path. While our Fermi liquid approach is not well con- 
trolled at high interaction, it has been proven to work at 
relatively large Vg which substantiates its use in the high 
density, large conductance regime kpl > 1. We therefore 
mostly show our theoretical results in the kpl > 1 metal- 
lic regime except for the lowest density (dashed) curve in 
Fig. 3(a). We emphasize that by adjusting the charged 
impurity density and by including disorder effects in the 
screening function, it is possible to get reasonable quan- 
titative agreement between our Boltzmann-Fermi liquid 
theory [Fig. 3(a)] and the experimental data (Fig. 2) 
even in the insulating low-density kpl < 1 regime (20j . 
We show our calculated results in Fig. [3] for various hole 
densities. In our calculation the overall resistivity scale 
depends on the unknown charge impurity density, but the 
qualitative trends in p{T,p) arise from basic aspects of 
the underlying scattering mechanism. Our calculated re- 
sistivity. Fig. [3ja), shows a non-monotonicity (as seen in 
Fig. [2] of the experimental data) where it increases with 
T at first and then decreases with increasing T after going 
through a maximum at a density-dependent temperature 
Tfn where dp/dT changes sign. This non-monotonicity is 
a non-asymptotic finite temperature "quantum-classical" 
crossover phenomenon, and arises from the competition 
between two contributions to transport [1], i.e. the tem- 
perature dependence of the dielectric screening function 
and the thermal averaging of the scattering rate. The 
screening effect decreasing with increasing temperature 
gives rise to increasing effective disorder with increasing 
T , and hence increasing p{T) with T. Based on the Boltz- 
mann theory, however, in the higher temperature classi- 
cal non-degenerate regime the thermal averaging of the 
scattering rate becomes quantitatively more important 
and shows stronger temperature dependence than screen- 
ing, which gives rise eventually to a p{T) decreasing with 
T for r > Tm- Thus, if the system is homogeneous and 
not strongly localized this crossover behavior appears in 
the transport measurements as can be seen in Fig. [51 
We note that as the density is reduced significantly be- 
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FIG. 3: (a) Calculated resistivities for p-GaAs quantum well 
system (a = 200 A) for various hole densities p =0.4, 0.5, 
0.6, 0.7, 0.9, 1.1 xlO^° cm"^ (from the top) with linearly 
screened charge impurity scattering. (The dashed curve is in 
the kpl < 1 regime.) (b) Calculated mobility as a function of 
hole density at r=50 mK. The dashed straight line indicates 



the /i ^ 
regime. 



p behavior. We only show mobility in the kpl > 1 
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low p = 9.0 X 10^ cm^^, kpl ~ 1 and direct comparison 
between theory and cxpcrmcnt becomes more problem- 
atic. Nevertheless, we point out that the theory repro- 
duces well all the qualitative features of the data in the 
metallic regime shown in Fig. [21 (1) An approximately 
(and upto) a factor of 2 increase in low-temperature p{T) 
with increasing T; (2) a maximum in p(T) at a density 
dependent Tm{p)', (3) a "high-temperature" decrease of 
p{T) ~ for T > Tm. One important difference be- 
tween the 2D n-GaAs electron system [1, 113, 1I| and our 
2D p-GaAs hole system is the very strong metallic tem- 
perature dependence (~ upto 200%) we observe (cf. Fig. 
2) in our hole sample versus the rather weak (~ 20% 
only) temperature dependence in the metallic phase of 

This strong (weak) 



11 



the n-GaAs electron system [IC 
metallic temperature dependence in the 2D hole (elec- 
tron) GaAs system is easily explained by the screening 
theory as arising from the much larger (by a fac- 
tor four) dimensionless screening parameter qxpl'^-kp in 
the hole system compared with the electron system due 
to the larger hole effective mass. In Fig. El^b) we show 
the calculated mobility as a function of hole density at 



T = 50 mK. At high densities (p > 1.0 x 10 



10 



the 



exponent a in /i oc is about 0.7, which is consistent 
with experiment (cf. Fig. |4|) indicating that the screened 
charged impurity is the main scattering source. 

In order to gain further insight into the transition to in- 
sulating behavior at low density we now focus on the den- 
sity dependence of the conduction properties at fixed low 
temperature. Figure 4(a) shows the dependence of mo- 
bility on density at fixed temperature T=47 mK. While 
the data for all densities cannot be described by a single 
power law behavior, at densities above ^2x10^" cm~^ 
the mobility approaches an approximate power law be- 
havior ^ ~ p^^'^ as seen in our theory. At low densities, 
p< 1 X 10^° cm~^, the power law changes to a much 
stronger function of density, pointing to a change in the 



FIG. 4: (a) Mobility vs. density at fixed temperature r=47 
mK. (b) Conductivity vs. density (solid squares) along with 
the fit generated assuming a percolation transition. The 
dashed line in (a) indicates the ji ~ p"'^ behavior. 
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FIG. 5: a) Variation of the derived critical density pc as a 
function of temperature, b) Variation of 5 with temperature. 



screening properties of the system at low carrier density. 
We emphasize that the quantitative and qualitative dis- 
crepencies between the lowest density data of Fig. 4a and 
the calculation of Fig. 3b indicate that the screening the- 
ory which captures the essential experimental features at 
higher densities does not fully describe the transition to 
insulating behavior. 

For high mobility n-type GaAs, Das Sarma et al. [llj 
interpreted the density driven crossover from metallic 
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conductivity to an insulating regime as a percolation- 
type transition. The importance of interactions as pa- 
rameterized by Tg for the MIT in GaAs remains an open 
question. In high mobility n-type 2D layers, the regime 
of r^'^-'lO can be reached at low density. Due to the higher 
effective mass in our 2DIIS, reaches approximately 45 
at low density in this study. It is an interesting ques- 
tion if a percolation type analysis can also be applied to 
our ultra high mobility 2DHS in the regime of large rg. 
We mention that there have also been various percola- 
tion anal yse s of 2D MIT in other n- and p-GaAs sys- 
tems (l^. l2l| for low mobility (a factor of 10 lower than 
our mobility) systems where the temperature dependent 
metallic behavior is much weaker than what we sec in 
our high-mobility system. 

To describe the transition of the resistivity from metal- 
lic behavior to an apparent insulating behavior at low- 
est densities, we fit our data using a percolation tran- 
sition model [i^l. The physical motivation for a per- 
colation transition is readily apparent. As the density 
of the 2DIIS is reduced, the system can no longer fully 
screen the random disorder potential. The 2DIIS breaks 
up into "puddles" of carriers isolated from one another 
by peaks in the disorder potential. Transport then takes 
place through percolating conduction paths in the dis- 
order potential network. In the percolation model, the 
conductivity near the percolation threshold density is de- 
scribed by: 

a{p) ^ A{p - p,)' 

where cr(p) is the density dependent conductivity, pc is 
the critical density, and 6 is the critical exponent de- 
scibing the transition. A is a constant of proportion- 
ality. The 2T> percolation exponent is expected to be 
1.31 ll|, |2lJ, I24I. The results of such an analysis ap- 
plied to our sample is shown in Figure 4(b). The ex- 
tracted exponent for our sample is (5=1.7±0.1 and the 
critical density Pc=0.38xl0^° cm-^iO.OlxlO^" cm'^ 
Formally, the preceding analysis is valid for (jp—pc) Pc- 
While the quality of the fit does not depend sensitively 
on the number of data points included, clearly the re- 
striction (p — Pc) <C Pc rnust be relaxed somewhat to 
make connection with the experimental data. Interest- 
ingly, the extracted conductivity parameters from the 
percolation analysis are consistent with known behav- 
ior of very high mobility n-type GaAs 2D systems as 
reported in ref. [ll] |. In particular, our extracted perco- 
lation exponent of 1.7 at r=47 mK is almost identical 
with the corresponding exponent value at a somewhat 
higher temperature 100 — 200 mK) in ref. ll|, which 
is consistent with equivalent values of the dimensionless 
temperature T/Tp for the two systems, where Tp is the 
Fermi temperature, and our critical density of 3 — 4 x 10^ 



cm ^ is consistent with the corresponding critical elec- 



tron density in ref. ll|. The temperature dependences 
of the fit parameters are detailed in figure 5. To con- 
clude, we find that the temperature dependence of the 
metallic phase is well-described by assuming the conduc- 
tivity to be limited by temperature and density depen- 
dent screened charged impurity scattering whereas the 
apparent low-density {r^ 3.8 x 10^ cm~^) transition from 
the metal to the insulating phase is well-described by a 
density-inhomogeneity driven percolation transition. 
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